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In April of 1992 Dr. Robert Stein visited Dr. D. Hathaway 
(SSL/MSFC/NASA) and the P.I. to discuss magnetohydrodynamic waves through 
the solar atmosphere. 

Dr. Marta Rovira, from the Institute of Astronomy at Buenos Aires, Argentina 
visited SSL/MSFC/NASA and the University for a two week period in April 1992 
to discuss solar physics and magnetic field research. 

September 1992 the P.I. had extensive discussions with the 
SSL/MSFC/NASA personnel (Dr. J. Davis) concerning the plans for the 
instrumentation development on board rocket projects. In particular the 
development of vector solar magnetograph on balloon flights. 

Dr. Robert Rosner was invited and accepted an invitation to returned to 
Huntsville in March 1993 to discuss further the present magnetohydrodynamic wave 
generation on the solar surface and its effects on the dynamics in the solar and 
stellar atmosphere. 

Dr. S. T. Wu and Dr. J. Davis again discussed the development of vector 
solar magnetograph on balloon flights as part of the instrumentation planning of 
rocket projects in June 1993. 

In October 1993 Dr. Chin-Chun Wu performed asstudy of slow shock 
evoluiton by numerical simulation modeling and gave a report to the P.I. These 
results are scheduled to be presented at the 1993 Annual Fall Meeting of the 
Amerian Geophysical Union. Some of the key results are included in the appendix. 



APPENDIX 


Numerical Study on the MHD Slow Shock Generation and 
Propagation in the Solar Wind 

S. T. Wu, C. C. Wu, (Center for Space Plasma and Aeronomic 
Research and Department of Mechanical Engineering, The 
University of Alabama in Huntsville, Huntsville, AL 35899; 
205-895-6413; e- mail wu@cspara.dnet.nasa.gov) 

S. T. Suess, Space Science Laboratory, NASA/MSFC, Huntsville, 

AL 35812 

J. K. Chao, Institute of Space Science, National Central University, 

Chung li, Taiwan, Republic of China 

It is well-known that most MHD shocks observed at 1 AU are MHD fast 
shocks. Only very limited number of MHD slow shocks are observed at 1 
AU. In order to understand why there are only a few MHD shocks 
observed at 1 AU we have performed a numerical study using an adaptive 
grid, unsteady, two-dimensional MHD model (Pantichob, Wu and Suess, 
1987 AIAA Paper 87-1218, Washington, DC) to investigate the MHD 
slow shock generation and propagation in the solar wind. In these 
numerical experiments a total of twenty-two cases of numerical 
calculations with various boundary perturbation are performed. These 
numerical results are summarized as follows: 

1. The forward slow shock (FSS) and reversed slow shock (RSS) can 
pass through each other and keep its own characteristics 

2. The second FSS will catch the first FSS and emerge into a stronger 
FSS. 

3. The FSS always disappears within a distance of 140 Rj (solar radii) 
from the sun when there is a forward fast shock propagating in front of 
them. 

4. When a FSS propagates behind a FFS, it shows the decreasing mach 
number of FSS. 

5. When a FSS propagates in front of a forward fast shock (FFS), it 
always will be caught by the FFS and destroyed by it. 

6. In all the tests we hve performed we have not discovered that the FSS 
(RSS) evolves into a FFS. 

This work supported by NASA Headquarters Grant NAGW-9. 
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NUMERICAL EXPERIMENTS FOR SIMULATION OF SHOCKS 


Q 

D 

H 

CO 

< 

O 

55 

w 


£ 

o 


c 

3 

cd 


cd 

x 

£ 

ca 

73 

C 

3 

u 

ca 


ca 

CN 

u- 

'i—* 

3 


CO 

CD 


> •« 


> 

H 

CL 

CD 

co 

u- 

cd 

CD 

£ 

cd 

cd 

CL 

X 

co 

CD 

u. 

cd 

cd 

u- 

CD 


C 

o 
+- * 
cd 

=3 


C 

CD 

co 

CD 

i— 

CL 

OD 


73 

od 

d 

CO 

£ 

00 

cd 

X 

4— » 

oo 

3 

E 

o 

II 

tcq 

> 

CD 

CD 

.3 

GO 

GO 

JJ 

CJ 

O 

X 

oo 

Q 

X 


cd 

cd 

u- 

CD 

C 

cd 

M 


73 

cd 

JO 

u- 

3 

ts 

CD 

CL 

cd 

JO 


IE 

cd 

G 

cd 

cd 

r; 

OX) 

cd 

E 


73 

C 

cd 


CL 

£ 


o 

CO 


00 

C 

CD 

E 

•5 

i 

CD 

C 

o 

0) 


c 

3 

CO 

C 

O 

a 


o 

CO 

a 

£ 

CD 

<D 


3 CO 

cd 3 


o 

Q- 

E 

o 

a 

-©- 

cd 


c 

o 


03 CD 

=§ € 

Jf c 

C cd 

“ X 

aj u* 

CO CD 

cd — 

CD Cd 

c E 

O oo 


CD 

73 

33 

‘E 

00 

cd 

E 

O 

co 

CD 

73 

u- 

O 

O 

£ 

CD 

i— 

cd 


id 

JO 


CD 
u- 

O 
C 

£> .2? 
a _>> 

CD 

> 


CL 

E 

CO 

CD 

> 


JO 

H 

CO 

» 

C 

CD 

C 

o 

CL 

E 

o 

o 


cd 

JO 

GO 

% 

ajT 

i— 

£ 

(D 

U* 

CD 

H 

CO 

C 

O 

cd 

JO 

u« 

3 

ta 

<D 

tO- 

co 

w 

c 

CD 

3 

o 

CL 

E 

0 

CD 

1 

-©- 

CD 

J= 


o 

'4— » 

c 

o 


c 

o 

3 

JO 

3 

tS 

CD 

CL 

CD 

JO 

■— * 

c 

o 

73 


3 

O 

co 

3 

CD 

*2 


— j= 


CO 

CO 

cd 

X) 

CD 

X 

-#— » 

c 

O 


73 

C 

cd 

rv 

u. 

> 

CL 

co' 

CD 

CD 

£ 

3 

i— . 

cd 

CL 

*cd 

a 

co 

>> 


CO 

3 

# o 

X 

cd 

> 


o •B 


cd •- 
X £ 

l— 

3 CO 

C 

o 


t: 

CD 

CL 

CD 

X 

CO 

CO 

cd 

CD 

CD 

> 


3 

C 

CD 

CD 

CO 

CO 


3 

X 

u- 

3 

CD 

CL 

c^ 

O 

CO 

P 


CD 

CO 

CD 

X 

W) 

c 

co 

D 


CO 

— t 

3 

CD 


CD 

CL 

X 

CD 

<N 

CN 

73 

CD 

E 


CD 

CL 

CD 

> 

3 

X 

CD 

£ 

o 

H 


CD > 


o 

> 


3 

.2? 

IE 

c 


o 

X 

co 

CO 

3 

CO 

CD 

x 

O 
0 JO 
3 

3 

CD 

X 

co 

W) 

c 

s 

CD 


> X 


J3 CL <4-^ 


o' 

CD 

CL 

X 


3 

Q. 

H 



co 

c 



CD 

73 

CD 

CD 

E 

U 

3 

‘a. 

3 

CL 

CD 

73 


CD 

3 

CD 

Ui 

3 

CO 

CO 
'*— » 


C 

X 

CD 

CL 

£ 

U-, 


o 

X 


CD 

co 

r- 

CL 

X 

3 

CD 

3 

CM 

C 

<N 

X 

CD 

CO 

p 

CD 

5 

X 

CD 

H 



\ 

Q. 


0.0 C J 

0 5 10 15 20 

Time( hours) 



o 10 20 30 40 

Time( hours) 


Figure 1. The Types of Perturbations. 






Table II Listing of 22 Cases of Six Different Types of 

Perturbations 



first perturbation 

second perturbation 

1 







Perturbation i 

Case 

4- 

P 0 

v r 

Vr 0 

V 

To 

tt 

£— 

P 0 

a 

i — 

r. 

Type 

i 

0.4 






I 

2 

0.8 






I 

3 


0.4 





I 

4 


0.8 





I 

5 



0.4 




I 

6 



0.8 




I 

7 

1.6 






II 

8 

3.2 






II 

9 


1.6 





II 

10 


2.0 





II 

11 



1.6 




II 

12 



3.2 




II 

mm 

0.4 

0.4 






IS 

0.4 


0.4 





15 


0.4 

0.4 





16 

5.0 

1.6 





IV 

17 

2.0 


5.0 




IV 

18 


1.6 

5.0 




IV 

19 | | 

0.4 | 



3.2 | 

V 

20 


0.4 


5.0 

1.6 j 

VI 

21 


0.4 

0.4 

5.0 

1.6 | 

VI 

22 


0.4 

0.4 ! 

1 

1.6 1 5.0 

VI 


where subscript "o" denotes the initial values at the reference 
point, which is at 28 R s ; superscript denotes the variable 
value of first perturbation; """ denotes the variable value of 
second perturbation at the low boundary which is at 28 R s for 
this study; 0.4 (0.8) means the perturbation decreasing to 40 
percent (80 percent) of its initial value; 1.6 (2.0, 3.2, 5.0) 
denotes the variable of the perturbation increasing to 1 60 (200, 
320, 500) percent of its initial value. 



















NUMERICAL RESULTS 
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Table III. Summary Results of Type I Perturbations 



















Table IV. Summary Resuits of Type II Perturbations 
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Pi 
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Pi 
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•i O 
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fi = 3 ' 2 
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1 FSS I 
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case 
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II 
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1.42 
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Pi 
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1.66 

11 

h - L6 
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1.47 
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at 
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10 
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Table V. 



Case perturbation 

13 £ = 0.4, £=0.4 
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Table VT. Summary Results of Type IV Perturbations 
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Figure 2. The evolution of radial velocity for case 1 , — = 0.4 
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Figure 3. The evolution of radial velocity for case 9, — = 1.6 
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II. Multiple Pulses (Square Waves) Perturbations (i.e. Type V and VI): 
Shock -Shock Interactions 
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II. Multiple Pulses (Square Waves) Perturbations (i.e. Type V and VI): Shock -Shock Interactions (Cont.) 
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Table VTI: Summary Results of Type V Perturbations 



I i FFS(B) 




I formed at 


i Mf, at 

second 

Case 

perturbation 

Time 

r 

ave. 

i 140 

180 

perturbation 



(hrs) 

(*.) 

Mf 

1 R, 

R. 

starts at (hrs) 

19. & 

11-0 4.11-32 

T. T. ~ 

49 

100.0 

1.36 

1.41 

1.55 

15 

19. b 

£ = 0.4; £ = 3.2 

69 

120.9 

1.48 

1.44 

1.51 

25 

19. c 

£ = 0.4; £ = 3.2 

81 

132.0 

1.50 

1.26 

1.52 

33 



FSS{ A+B) 




formed at 


M 

, at 

second 

Case 

i perturbation 

Time 

i r 

ave. 

140 

180 

perturbation 



(hrs) 

(R.) 

M f 

R. 

R. 

starts at( hrs) 

19.a 

£ = 0.4; £ = 3.2 

49 

93.5 

2.23 

2.15 

2.24 

15 



FSS(A) 




formed at 


Af« at 

second 

Case 

perturbation 

Time 

r 

ave. 

140 

180 

perturbation 



(hrs) 


Mf 

R, 

R, 

starts at(hrs) 

19. a 

£ = 0.4; £ = 3.2 

13 

33.1 

1.36 

29 

59.0 

15 

19.b 

£ = 0.4; £=3.2 

13 

32.7 

1.623 

81 

145.9 

25 

19.C 

£ = 0.4; £ = 3.2 

13 

33.1 

1.72 
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184.5 
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1 1 FSS(B) I I 



formed at 1 


Mf, at 1 

second 

Case 

perturbation 

Time 

r 

ave. 

140 

180 

perturbation 



(hrs) 

(*,) 

Mf 

R. 

R, 

starts at (hrs) 

19. b 

£ = 0.4; £=3.2 

29 

35.3 

1.80 

69 

114.8 

25 

19. c 

£ = 0.4; £ = 3.2 

37 

35.1 

1.69 

113 

184.5 

33 


Where (A) denotes the shock generated by the first perturbation; (B) denotes 
the shock generated by the second perturbation; (A+B) denotes the merged 
shock from two shocks which were generated by the first and second 
perturbations. 
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200 
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R f 

R„ 
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25 
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3.14 
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1.6 
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48.6 

3.46 

4.56 

2.35 

= 5.0 
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48.1 

3.16 

4.44 

2.47 


FFS(B) 


formed 


M, 

at 
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60 
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25 
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= 5.0 

25 
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Table Vmb. Summary Results c 
Starts at t = 33 hrs. 
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VI Perturbation; Second Perturbation 


j FFS(B) 

I formed at I M, at 



Time 

(hrs) 
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(*.) 

ave. 

M. 

60 

R. 

200 

R. 

= 1.6 

37 

37.4 

1.81 

1.38 

2.30 

= 1.6 

49 

66.3 

1.83 

1.10 

1.95 

= 5.0 
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1.20 

2.58 
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FSS(B) 


1 formed 


disappeared 


Time 
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ave. 

Time 
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(hrs) 

(R.) 

Af, 

(hrs) 

(*.) 


1.6 41 45.23 1.97 85 147. 

1.6 37 35.86 1.45 73 116.8 

5.0 37 36.33 1.53 73 119.2 

FSS(A) 

formed disappeared 

Time r ave. Time r 
(hrs) (.R.) Af. (hrs) (fl.) 

1.6 13 31.1 2.71 85 154.8 

1.6 13 31.9 2.58 101 187.1 

5.0 11 29.6 2.00 89 162. 
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Figure 4. The evolution of radial velocity for case 19a, -£-=0.4, £—=3.2 
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Figure 5. The evolution of radial velocity for case 19b, -£-=0.4, £-=3.2 
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Figure 6. The evolution of radial velocity for case 19c, 


-£ 1 = 0 . 4 , £ 1 = 3.2 
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Figure 7 . The evolution of radial velocity for case 21b, 
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Figure 8. The evolution of radial 


velocity for case 21a, = 0.4, 
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On the basis of these twenty-two numerical experiments, we found: 
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